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ABSTRACT 

One of the central goals of the Spitzer Legacy Project "From Molecular Cores 
to Planet-forming Disks" (c2d) is to determine the frequency of remnant circum- 
stellar disks around weak-line T Tauri stars (wTTs) and to study the properties 
and evolutionary status of these disks. Here we present a census of disks for a 
sample of over 230 spectroscopically identified wTTs located in the c2d IRAC 
(3.6, 4.5, 4.8, and 8.0 /xm) and MIPS (24 /im) maps of the Ophiuchus, Lupus, 
and Perseus Molecular Clouds. We find that ~20% of the wTTs in a magnitude 
limited subsample have noticeable IR-excesses at IRAC wavelengths indicating 
the presence of a circumstellar disk. The disk frequencies we find in these 3 re- 
gions are ~3-6 times larger than that recently found for a sample of 83 relatively 
isolated wTTs located, for the most part, outside the highest extinction regions 
covered by the c2d IRAC and MIPS maps. This discrepancy in the disk fraction 
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of these two different groups of wTTs (on cloud vs. off-cloud targets) supports 
the idea that samples of wTTs distributed around molecular clouds (nominally 
1-10 Myrs old) represent a somewhat older population of stars. The disk frac- 
tions we find are more consistent with those obtained in recent Spitzer studies 
of wTTs in young clusters such as IC 348 and Tr 37. From their location in the 
H-R diagram, we find that, in our sample, the wTTs with excesses are among 
the younger part of the age distribution. Still, up to ~50% of the apparently 
youngest stars in the sample show no evidence of IR excess, suggesting that the 
circumstellar disks of a sizable fraction of pre-main-sequence stars dissipate in 
a timescale of ~1 Myr. We also find that none of the stars in our sample ap- 
parently older than ~10 Myrs have detectable circumstellar disks at wavelengths 
< 24 fim. Our result on the survival time of primordial disks is very similar 
to those obtained by studies based on K-band identified inner disks (r < 0.1 
AU). Also, we find that the wTTs disks in our sample exhibit a wide range of 
properties (SED morphology, inner radius, Ldisk/^*, etc) which bridge the gaps 
observed between the cTTs and the debris disk regimes. This strongly suggests 
that wTTs disks are in fact the link between the massive primordial disks found 
around cTTs and the debris disks observed around main-sequence stars. 

Subject headings: infrared: stars — planetary systems: protoplanetary disks - 
stars: pre-main sequence — open clusters and associations: individual (IC 348) 

1. Introduction 

Over the last couple of decades, it has been clearly established that circumstellar disks 
are an integral part of the star formation process. Even though there is currently no direct 
evidence that planets actually grow from circumstellar material, it has become increasingly 
clear that disks are potential birthplaces of planets because their masses, sizes, and com- 
positions are consistent with those of the assumed pre-planetary solar nebulae (Hillenbrand 
2003). More recently, the discovery of exo-planets orbiting nearby main sequence stars has 
confirmed that the formation of planets is a common process and not a rare phenomenon 
exclusive to the Solar System. 

Direct detection of forming planets is well beyond our current capabilities and observing 
molecular hydrogen, which largely dominates the mass of primordial disks, is particularly 
challenging (Thi et al. 2001, Richter et al. 2002). However, the thermal emission from 
circumstellar dust is much easier to detect and study. For this reason, the study of the 
evolution of circumstellar dust has been a natural first step toward providing observational 
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constraints on planet formation theories. 

Strom et al. (1989) studied a sample of 83 classical T Tauri stars (cTTs) and weak-line 
T Tauri stars (wTTs) located in the Taurus-Auriga star-forming region in order to determine 
the fraction of objects with K-band (2.2 /im) and IRAS excesses indicating the presence of a 
circumstellar disk. WTTs are low-mass pre-main sequence stars that occupy the same region 
of the H R diagram as cTTs but do not show clear evidence of accretion. The distinction 
between the two is usually made based on the Ha equivalent width (EW). The H a EW 
of cTTs is > 10 A, while the Ha EW of wTTs is < 10 A. Since there is a very strong 
correlation between spectroscopic signatures of gas accretion and the presence of near-IR 
excess (Hartigan et al. 1995), most CTTS show near-IR excess while most wTTs lack such 
an excess. Strom et al. (1989) found that 60% of their stars younger than 3 Myrs showed a 
K-band excess, indicating the presence of a circumstellar disk, while only 10% of the stars 
older than 10 Myrs years did so. Based on these numbers, they estimated a disk dissipation 
timescale of < 3-10 Myrs and claimed that their result was, at the time, "the best astrophys- 
ical constraint on the time available for planet building" . It has been argued that individual 
star-forming regions such as Taurus lack the intrinsic age spread necessary to investigate 
the dissipation timescale of circumstellar disks from individually derived ages (Hartmann 
2001). However, similar disk lifetimes studies based on the disk frequency of clusters with 
different mean ages and extending to the 3.4 fxm L-band (e.g. Haisch et al. 2001) have lead 
to results similar to those presented by Strom and collaborators (see Hillenbrand (2006) for a 
review on the frequency of near-IR excesses based on a sample of ~3000 PMS). The K-band 
excess, when used as a disk indicator, is only sensitive to dust in the innermost part of the 
disk; therefore, K-band studies only constrain the dissipation timescale of a region of the 
disks that is much closer to the star than the locations corresponding to the orbits of any 
of the planets in the Solar System (Mercury's semi-major axis is 0.39 AU). The dissipation 
timescale of the dust in the planet-forming regions might or might not be the same. Since 
the dynamical timescale is shorter and the surface density is higher closer to the central 
star, circumstellar disks are expected to evolve from the inside out. Most wTTs, which by 
definition present little or no evidence for accretion, also show little or no near-IR excess 
(Hartigan et al. 1995). However, even after the inner accretion disk has dissipated, it is 
entirely possible that wTTs to still have enough material at larger radii to form terrestrial 
and giant planets. In fact, millimeter-wave observations show that at least 10 % of wTTs 
have disks with estimated masses in the 10 _1 -10 -3 M range (Osterloh & Beckwith 1995, 
Andrews and Williams 2005). 

While the existence of planets around a significant fraction of all MS stars has been 
verified (e.g., Marcy & Butler 1998), the fundamentals of the planet formation process still 
remain open questions, especially for giant planets. There are currently three main theories 
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for the formation of giant planets: core accretion, gravitational instability, and hybrid models 
that combine aspects of both theories. See Lissauer & Stevenson (2007) and Durisen et al. 
(2007) for two recent reviews of the core accretion and gravitational instability models and 
a discussion of the many upstanding questions. 

Although it is unlikely that an observational estimate of the disk's dissipation timescale 
by itself can distinguish between the competing theoretical models mentioned above, esti- 
mates of the dissipation timescale of the planet forming regions can impose valuable con- 
straints on current theoretical models. In order to probe the planet-forming regions of disks 
around pre-main sequence (PMS) stars, mid- and far-IR observations are required. Unfortu- 
nately, these spectral regions are not easily observable from the ground, and past space IR 
telescopes such as IRAS and ISO were only sensitive enough to detect very bright optically 
thick disks in low-mass stars at the distance of nearest star-forming regions. These instru- 
ments lacked the sensitivity needed to detect the modest IR-excesses expected for optically 
thin disks and faint optically thick disks. Thus, Spitzer's sensitivity is required to establish 
whether most wTTs have optically thin disks, disks with inner holes, which are too cold to 
be detected in the near-IR from the ground, disks too faint to be detected by IRAS and ISO, 
or no disks at all. One of the main goals of the Spitzer Legacy Project "From Molecular 
Cores to Planet-forming Disks" (c2d; Evans et al. 2003) is to determine whether or not 
most wTTs have circumstellar disks and to characterize their properties and evolutionary 
status. Preliminary results from the c2d Legacy Project (Padgett et al. 2006, P06 hereafter) 
showed that disks are rare (~6%) among the population of wTTs distributed around nearby 
molecular clouds. However, other recent Spitzer studies have reported significantly larger 
disk fractions (~30%) among wTTs in young clusters such as IC 348 and Tr 37 (Lada et 
al. 2006; Sicilia-Aguilar et al 2006). Here we study a sample of over 230 spectroscopically 
identified wTTs located in the c2d IRAC and MIPS maps of the Lupus, Ophiuchus, and 
Perseus Molecular Clouds in order to investigate the frequency of circumstellar disks as a 
function of stellar age. In Section 2, we describe the c2d survey of molecular clouds and 
our sample of wTTs. In Section 3, we identify IR excesses and investigate the properties 
of their disks. In Section 4, we compare our results to previous Spitzer results and discuss 
the evolutionary status of wTTs disks. Also in Section 4, we derive the ages of the wTTs in 
our sample from their location in the H-R diagram. We investigate the disk frequency as a 
function of stellar age and use our results to impose constraints on the timescale for planet 
building. Finally, our conclusions are summarized in Section 5. 



- 5 - 



2. Observations 

2.1. C2D large molecular clouds and GTO observations 

As part of the c2d Legacy Project, Spitzer has mapped 13.0 sq. deg. of 3 nearby 
star-forming regions, Perseus, Ophiuchus, and Lupus, with IRAC (3.6, 4.5, 5.8, and 8.0 /im) 
and 22.1 sq. deg. with MIPS (24 /xm). MIPS 70 and 160 fim observations were also taken, 
but due to sensitivity considerations, we do not include these observations in most of our 
analysis. The IRAC maps consist of two dithers of 10.4 sec observations, each obtained 
at two epochs (41.6 sec total) separated by several hours. The second epoch observations 
were taken in the High Dynamic Range mode, which includes 0.4 sec observations before 
the 10.4 sec exposures, allowing photometry of both bright and faint stars at the same time. 
MIPS observations were taken with the fast scan mode, also in two different epochs of 15 
sec exposures each. See Jorgensen et al. (2006) and Young et al. (2005) for a detailed 
description of the observing strategy used for the c2d IRAC and MIPS survey of nearby 
molecular clouds. 

In addition to the data from the c2d legacy project, we use observations of IC 348 taken 
as part of the IRAC and MIPS Guaranteed Time Observer (GTO) programs (Program ID 
36 and 58, respectively). The IRAC GTO observations cover a 15'x 15' field of view centered 
in the cluster and consists of two pairs of 8 dithers of 96.8 sec exposures for the 3.6, 4.5, and 
5.8/zm observations (e.g. 1549 sec exposures per pixel). The 8.0/zm observations consist of 
four pairs of 8 dithers of 46.8 sec exposures 0. The MIPS 24 /im GTO observations of IC 348 
were taken in the medium scan mode resulting in an average exposure time of 80 seconds per 
pixel. For consistency, we processed the Basic Calibrated Data from the GTO programs and 
produced point source catalogs using the c2d pipeline (Evans et al. 2005). The c2d pipeline 
uses the the c2d mosaicking/source extraction software c2dphot (Harvey et al. 2004), which 
is based on the mosaicking program APEX developed by the Spitzer Science Center, and 
the source extractor Dophot (Schechter et al. 1993). 

Flux uncertainties in c2dphot are calculated in a standard way from a numerical esti- 
mate of the Hessian matrix (Press et al. 1992; Silvia 1996). This procedure for estimating 
uncertainties, although statistically correct, appears to underestimate the uncertainty as 
measured by the repeatability of flux measurements of the same objects at different epochs. 
For bright sources, there appears to be a random error floor to the best uncertainty possible 
with our observing techniques of 0.05 mag for the IRAC bands and 0.09 mag for the MIPS 



1 The longest integration time of the 8.0 /im array, nominally 100 sees, consists of 2 exposures of 46.8 sec 
each (See IRAC handbook, http://ssc.spitzer.caltech.edu/irac/dh/). 
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bands. The absolute calibration uncertainties are not included in our uncertainties. They 
are 5 and 10% for IRAC and MIPS, respectively (see data handbook for the instruments). 
As in most Spitzer surveys, the intrinsic sensitivity of the c2d observations is not uniform 
across the clouds due to variations in the total exposure time at different positions in the 
sky, in the amount of extended cloud emission, and in the source confusion level. Based on 
the cumulative fraction of all sources detected in both epochs of the c2d observations, the 
overall 90 % completeness levels of the c2d survey have been estimated to be 0.07, 0.12, 0.5, 
0.4, and 1.0 mJy for 3.6, 4.5, 5.8, 8.0, and 24 /jm. See Evans et al. (2005) for a detailed 
discussion of the uncertainties and sensitivity limits of the c2d survey. 

The IRAC observations are sensitive enough to allow robust detections of stellar pho- 
tospheres in all four IRAC bands for our entire sample of wTTs in Ophiuchus and Lupus 
(distance ~125 pc and ~150-200 pc, respectively) and in ~85% of the objects in Perseus 
(distance ~320 pc). MIPS 24 /zm observations are not deep enough to reach the stellar pho- 
tosphere of some low mass objects (especially in Perseus, which is the most distant cloud we 
consider), but, in general, they are deep enough to detect optically thick disks in our entire 
sample. The different completeness levels of our disk census due to sensitivity considerations 
are discussed in more detail in Section 13.1.11 

2.2. Sample Selection 

The c2d maps contain several hundred young stars identified by their X-ray and Ha 
brightness that have been spectroscopically classified as wTTs stars. Our sample was se- 
lected from the literature and is distributed as follows: 69 objects in Ophiuchus (Bouvier 
& Appenzeller 1992; Martin et al. 1998), 33 in Lupus (Hughes et al. 1994; Krautter et al. 
1997), and 130 objects in Perseus (Luhman et al. 2003). The Lupus and Ophiuchus objects 
are distributed across the cloud maps, while the targets in Perseus are located exclusively 
in the IC 348 cluster. Of the 33 Lupus objects, 27 are located in the c2d maps of Lupus III 
and 6 lie within the c2d maps of Lupus I. 

Spitzer SEDs of the wTTs in IC 348 have already been presented by Lada et al. (2006, 
L06 hereafter). We include these objects in our sample because they increase the statisti- 
cal significance of our results and allow us to compare a clustered population to the more 
distributed population of stars in Lupus and Ophiuchus. Also, we adopt a different disk 
identification criterion than L06, which results in a lower disk fraction than that obtained 
by the Lada group (see Section f4. 1.21) . Furthermore, L06 adopt a single age (2-3 Myrs) for 
the stars in the IC 348 cluster and do not attempt to study the disk fraction as function of 
age, which is one of the main goals of our paper. 
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All the objects in our sample, which are listed in Table 1, have known spectral types and 
small Ha EWs. The spectral types are necessary to estimate stellar ages from the position 
of the targets in the H-R diagram and the contribution of the stellar photospheres to the 
observed SEDs, while the Ha EWs are required to establish wTTs status. The nominal 
division between cTTs and wTTs is Ha EW = 10 A; however, since the Ha EW due to 
chromospheric activity alone can reach this value for late M stars (e.g., Martin 1997), we 
have included in our study 18 M2-M7 stars ( < 8 % of our sample) with Ha EW up to 15 
A. Also, ~8% of the stars in our sample show Ha in absorption rather than in emission. 

We note that even though there is a strong correlation between Ha emission and other 
accretion signatures such as optical veiling, a single-epoch low-resolution measurement of 
Ha equivalent width is not enough to rule out active accretion for at least two reasons. 
First, even when a narrow range of spectral types is considered, the distribution of Ha EWs 
of T Tauri stars does not show a clear gap between accreting and non-accreting objects 
(e.g., the Ha EWs of weakly accreting PMS stars overlap with those of chromospherically 
active non-accreting stars). Second, accretion itself is a highly variable process, and some 
objects constantly move across the wTTs-cTTs Ha EW boundary. Therefore, we consider 
our sample to be composed of mostly non-accreting objects but do not rule out the presence 
of some actively accreting interlopers. 



2.3. Complementary data 

In order to construct more complete SEDs of our wTTs, we have collected the 2MASS 
photometry for our entire sample and the V, Rc, /p-band photometry reported by Hughes et 
al. (1994) and Wichmann et al. (1997) for 24 of our 33 Lupus objects. Also we have obtained 
our own VRqIc optical observations for 52 of our 69 Ophiuchus targets and Rc and Ic-band 
observations for 115 of the 130 objects in IC 348 with the 0.8 m telescope at McDonald 
Observatory. The Ophiuchus targets were observed in 7 different 46'2 x 46f2 fields of view 
during the photometric nights of June 20* A -21 S * 2005 with exposures times of 30, 50, and 100 
seconds for the V, Rc, and Ic-band respectively. The objects in IC 348 were observed in a 
single field of view with 200 and 150 second exposures in the Rc and Ic-band, respectively. 
In addition to the program stars, on each night, 3 fields of Landolt standards (~5 standards 
per field) were observed at different airmasses. The seeing ranged from 1.5" to 2" when the 
observations were made. We reduced the data and perform aperture photometry using the 
standard IRAF packages CCDRED and DAOPHOT. We used a 5.4" (4-pixel) aperture and 
a sky annulus with inner and outer radii of 16.2" and 22.95" respectively. 

The rms scatter of the photometric solutions applied to the programs stars was < 0.02 
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mag in all three filters. We adopt a conservative minimum photometric error of 0.03 mag. We 
report the magnitudes and the uncertainties for all the objects with estimated photometric 
error less than ~0.2 mag. The Spitzer photometry for our entire sample is listed in Table 1. 
The non- Spitzer data: optical and 2MASS photometry, along with the spectral types and 
H a equivalent widths from the literature, is listed in Table 2. 

3. Results 

3.1. Disk Identification 

In order to identify the stars with disks, we compare the extinction-corrected Spitzer col- 
ors of our targets, to those predicted by NextGen Models (Hauschildt et al. 1999), convolved 
with the Spitzer bandpasses, for the photospheres of stars of the corresponding spectral types. 
The broader the wavelength baseline of the color used, the larger is the expected excess of the 
stars with disks; therefore, the available color that provides the most clear disk identification 
is [3. 6] -[24]. However, since both Spitzer's sensitivity and the photospheric fluxes decrease 
with increasing wavelength, not all sources are detected at wavelengths longer than 5.8 /xm. 

In Figure la, we plot [3.6] vs EX([3.6]-[24]), for the 98 stars with 3.6 and 24 fim fluxes 
available, where EX([3.6]-[24]) is ([3.6]-[24]) bs -([3.6]-[24]) m «w, ([3.6]-[24]) bs are the 
observed colors corrected for extinction, and ([3.6]-[24])ModeJ are the photospheric colors 
predicted by the NextGen Models. We estimate the extinction, Ay, using Ay = 4.76 x 
E(R — I) = 4.76 x ((Rc — Ic)obs — (Rc — Ic)o)- The intrinsic stellar colors, (Rc-Ic)o, come 
from Kenyon & Hartmann (1995). For objects without Rc and lc fluxes available, we use 
Ay=5.88xE(J-Ks). We note that this will result in an overestimated extinction for objects 
with significant K-band excess. The extinction at 8 and 24 fxm are calculated according to 
Table 3. 

We consider objects with [3.6]-[24] < 0.7 to be those without excess whose emission 
arises solely from the stellar photosphere. The mean EX([3.6]-[24]) value for this group is 
not zero but 0.07 mag with a I- a dispersion of 0.17 mag. The 0.07 mag offset of the observed 
stellar photospheres with respect to the models is probably due to a combination of the 
systematic errors in the absolute flux calibrations, the stellar models, and the extinction 
corrections. In Figure la, we treat this offset by subtracting 0.07 mag from all the EX([3.6]- 
[24]) values. We find that 40 objects have a [3.6]-[24] excesses larger than 5-cr. These are 
very robust disk identifications. One object, RXJ1622. 6-2345, has [3.6]-[24] excess just over 
3-er. We consider this object to be a good disk candidate, but warn the reader of its lower 
significance. Since the c2d MIPS maps cover a larger area than the IRAC maps, 24 fim is 
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the only Spitzer flux available for 5 stars in our sample. In these cases, we use [K]-[24] colors 
for disk identification. None of the 5 objects for which MIPS 24 /j,m is the only available 
Spitzer flux show a significant K-24 /xm excess. Of the 127 stars without measured [3. 6] -[24] 
or [K]-[24] colors, 112 have [3.6]-[8.0] colors available. Figure lb is analogous to Figure la, 
but here we plot [3.6] vs EX([3.6]-[8.0]) for all the stars with measured [3.6]-[8.0] colors, 
including the ones from Figure la, which are shown as open diamonds. Following Cieza & 
Baliber (2006), objects with [3.6]-[8.0] < 0.7 are considered stellar photospheres. In this case, 
the mean color offset with respect to the models is 0.05 mag. The standard deviation of the 
stellar photospheres is 0.16 mag, but the error clearly increases with decreasing brightness. 
We find that only 6 objects not detected at 24 fxm show a clear evidence (> 5-cr) for 8 /im 
excess. Two objects, IC348-76 and IC348-67, show excesses between 3 and 5-cr. Since the 
SEDs of these two objects show a hint of IR-excess at 4.5 and 5.8 jum, we consider these two 
objects to be good disk candidates. We find that none of the 15 objects without [3.6]-[24], 
[Ks]-[24], or [3.6]-[8.0] colors available shows a significant [3.6]-[5.8] excess. 

The SEDs of the wTTs disks in Lupus, Ophiuchus, and IC 348 are shown in Figures 2, 
3, and 4, respectively. The open squares represent the observed optical, 2MASS, IRAC and 
MIPS-24 fim fluxes, while the dots correspond to the extinction corrected values. The Ay's 
are estimated as described above for Figure 1, while the extinctions at other wavelengths are 
calculated according to Table 3. NextGen model photospheres, corresponding to published 
spectral types and normalized to the extinction-corrected J-band, are shown for comparison. 

3.1.1. Disk Census Completeness 

Following the procedure described above, we identify a total of 46 wTTs disks and 3 
disk candidates. For definiteness, we consider these 3 disk candidates to be real in the rest 
of the paper. However, we note that the completeness of our disk census is lower for our IC 
348 sample than it is for the sample of Lupus and Ophiuchus wTTs. 

The IRAC observations are sensitive enough to allow robust detections of stellar pho- 
tospheres in all four IRAC bands for wTTs in the sample from Lupus and Ophiuchus (102 
objects in total). Also, the 5 objects in Lupus and Ophiuchus that fall outside the c2d IRAC 
maps all have 24 /im fluxes consistent with stellar photospheres, which rules out the presence 
of any significant IRAC excess in these objects. Therefore, the census of IRAC excesses in 
our sample of Lupus and Ophiuchus wTTs is likely to be complete (i.e. only 17 out of the 
102 wTTs in Lupus and Ophiuchus have significant IRAC excesses). At 24 /im, we detect 
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80% (82/102) of the Lupus and Ophiuchus objects. Thus, it is possible that some of the 
20 wTTs from the Lupus and Ophiuchus sample that are not detected at 24 jum have SEDs 
that start to diverge from their photospheres longward of 8.0 /im but remain bellow our 
sensitivity at 24 /im. We find that all the objects with IRAC excess are detected at 24 /im, 
and that only 2 of the 82 objects detected at 24 /im have excesses that only become evident 
at this wavelength. 

In our IC 348 sample, 11% (14/130) of the objects are not detected at 8.0 /im and 82 
% (107/130) are not detected at 24 /im. The longer integration times of the GTO IRAC 
observations of IC 348 with respect to the c2d observation of Lupus and Ophiuchus (1548 
vs. 46.8 sees) more than compensates the effect of the distance (320 vs. 125-200 pc) on 
the expected sensitivities. However, the high background of the IC 348 cluster becomes the 
limiting factor for the detection of faint sources at 8.0 /im. At 24 /im, the depth of the 
GTO observations only partially compensates for the greater distance of IC 348, and the 
effect of the background becomes even larger than it is at 8.0 /im. The combination of these 
two factors explains the very low detection rate of IC 348 members at 24 /im. None of the 
14 objects without a 8.0 or 24 /im detection shows a significant 5.8 /im excess that would 
indicate the presence of a circumstellar disks. However, some of these 14 wTTs could have 
SEDs that start to diverge from their photospheres longward of 5.8 /im but remain bellow 
our sensitivity at 8.0 and 24 /im. Finally, we find that 91% (21/23) of the objects detected 
at 24 /im also show significant IRAC excesses, which is consistent with the results for Lupus 
and Ophiuchus. Obviously, in IC 348, we are likely to have missed most of the stars for 
which the onset of the IR-excess occurs longward of 8.0 /im. 



3.1.2. Disk fraction statistics 

In section [3TT1 we tried to present a disk census as complete as possible given the available 
data. However, since the sensitivity of our disk survey is not uniform across all wavelengths 
and varies from region to region due to distance and background level effects, we derive disk 
fraction statistics from the number of objects with excess at IRAC wavelengths (IRAC disk 
fraction, hereafter) in a magnitude limited subsamples that is not likely to be affected by 
sensitivity variations. 

As discussed in the previous section, we consider our census of IRAC excesses in our 



2 by detect we mean that the flux of the object has been measured regardless of whether or not the star 
has an excess. All the objects have been observed at 24 /im;therefore, a non detection implies that the flux 
of the object is below our 24 /im sensitivity limit. 
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sample of Lupus and Ophiuchus wTTs to be complete. Excluding the disks that are only 
detected longward of 10 /iin (ROXs 36, and RX J1622. 6-2345), we derive the following IRAC 
disk fractions: 27%±7% (9/33) for Lupus, 13%±4% (9/69) for Ophiuchus. For IC 348, we 
restrict our sample to the 96 objects with 3.6 /im fluxes greater than 3.2 mJy, a level above 
which we also detect at 8.0 fim every object in our sample. We find that 22 of these 96 
objects show significant IRAC excess, which leads to a disk fraction of 23±4%. Combining 
the objects in Lupus, Ophuichus, and IC 348, we derive an overall IRAC disk fraction of 
20±3% (40/198). The IRAC disk fractions of wTTs in Lupus and Ophiuchus bracket that 
of the IC 348 cluster. This suggest that the disk fraction in IC 348 wTTs targets is not 
strongly affected by the cluster environment. 

At 24 fim, the sensitivity of our survey is significantly less uniform than at IRAC 
wavelengths. We find that even some of our brightest objects lack 24 /im detections due 
to the strong background emission surrounding them. This prevent us from deriving robust 
MIPS disk fractions statistics even from a magnitude limited subsample. 

3.2. Disk Properties 

3.2.1. Selected Color- Color Diagrams 

In this section, we place our sample of wTTs in Spitzer color-color diagrams as our first 
attempt to explore the diversity of wTTs disks. For comparison, we also include in these 
diagrams two samples of 83 and 66 cTTs from Taurus and IC 348, respectively. We used the 
Ha EW from Luhman et al. (2003) to select the sample of cTTs in IC 348. For these objects 
we use our own photometry. The Spitzer fluxes of the Taurus cTTs come from Hartmann 
et al. (2005). In Figures 5 and 6, we show an Ay=l0 extinction vector based on Table 
3; however, we note that for most of the stars in our sample the extinction is significantly 
smaller. The mean Ay we derive for our Lupus, Ophiuchus, and IC 348 wTTs samples are 
0.64, 2.7, and 1.7 mag respectively. Also, 93% of our sample has Ay < 5 mag, and the 
remaining 7 % has 10 mag < Ay < 5 mag. 

Figure 5a shows the [3.6]-[24] vs. [3.6]-[8.0] colors of our sample of wTTs stars. Based on 
the colors shown, we identify three different groups. The first group consists of objects with 
[3.6]-[24] < 0.7 and [3.6]-[8.0] < 0.4, which are consistent with bare stellar photospheres. 
The second group consists of stars with [3.6]-[24] > 0.7 and [3.6]-[8.0] < 0.4 (i.e. their 
SEDs start to diverge from their stellar photospheres longward of 8.0 /im). The 4 stars in 
this group are RXJ1622.6-2345, ROXs 36, IC348-56, and IC348-124. In Section KIM we 
find that the objects in this second group have the lowest fractional disk luminosities of 
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the sample. In fact, all four objects have Ldisk/L* < 10~ 3 , which suggest the presence of 
optically thin disks. The last group of objects has [3.6]-[24] > 0.7 and [3.6]-[24] > 0.4 (i.e. 
they show excess at both IRAC and MIPS wavelengths). Most of these objects are likely to 
be optically thick disks. Figure 5b shows our sample of wTTs combined with the sample of 
cTTs. We note that the cTTs populate exclusively the region of the diagram we associate 
with optically thick disks and that both wTTs and cTTs populations are very well mixed in 
this region of the diagram. 

In Figure 6a, we show the [3.6]-[4.5] vs. [5.8]-[8.0] colors of wTTs. In general, the 
populations of stars with and without an excess are clearly separated. Stars in the upper 
right corner of the diagram have both 4.5 and 8.0 /xm excesses. Stars in the lower right 
corner of the diagram are stars with 8.0 fim excess but no 4.5 /xm excess. These objects are 
usually interpreted as "transition disks" with inner holes (e.g. Calvet et al. 2002; D'Alessio 
et al. 2005). In fact, Sicilia-Aguilar et al. (2006) define transition disks as objects with 8.0 
fim excess but no 4.5 fim excess. They find that only ~10% of the T Tauri disks in the 
4-Myrs old cluster Tr 37 fall into this category. In contrast, we find that ~30 % of the wTTs 
are in fact "transition disks" according the working definition stated above. This is in good 
agreement with the idea that wTTs disks represent a more evolved evolutionary state than 
cTTs disks. In Figure 6b, we combine wTTs and cTTs in the same color-color diagram. 



3.2.2. EW(H a ) dependence? 

Since the distinction between cTTs and wTTs is based on the H a EW, and cTTs show 
a much higher disk fraction than wTTs, we investigate the dependence of the IRAC disk 
fraction on Ha EW within our sample of wTTs. We restrict our analysis to the 198 objects 
in the magnitude limited subsample discussed in Section 13.1.21 Figure 7a shows a histogram 
of the Ha EW of stars with and without a disk, while Figure 7b shows the disk fraction 
vs. Ha EW. We find that (1) the disk fraction is highly correlated with Ha, and (2) the 
disk fraction is a smooth function of Ha EW, increasing from ~5% for the stars with Ha 
observed in absorption to >50 % for the stars with 15 A > Ha EW > 10 A. 

Since a strong correlation between the presence of a disk and Ha emission of chromo- 
spheric origin would be difficult to explain, the correlation observed in Figure 7 is likely to 
be driven by the presence of weakly accreting stars in our sample of wTTs. Muli-epoch 
high-resolution spectroscopy of the sample of wTTs with disks would be highly desirable to 
establish, from their Ha velocity profiles, which objects are in fact actively accreting (White 
& Basri 2003). 
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3.2.3. Disk fraction dependence on spectral type? 

Lada et al. (2006) find that the "optically thick" disk fraction in the entire populations 
of PMS stars (cTTs/wTTs) in IC 348 seems to be a function of spectral type. They argue 
that the disk fraction peaks around K6-M2 stars, which at the age of the cluster corresponds 
to stars with masses similar to that of the Sun, and conclude that circumstellar disks might 
last longer around solar type stars than around both less and more massive stars. Similarly, 
Carpenter et al. (2006) studied a sample of over 200 PMS star in the 5 Myrs old Upper 
Scorpius OB association with masses between ~0.1 and 20 M Q and found that stars with 
K-M spectral types have a significantly larger disk fraction than stars with G-B spectral 
types. 

In order to investigate the possibility of a disk fraction dependence on spectral type 
for wTTs, we restrict our analysis to the 198 objects in the magnitude limited subsample 
discussed in Section 13.1.21 and divide our sample into the same 3 spectral type bins studied 
by Lada et al. (2006). We find the following IRAC disk fractions for the 3 bins: 28%±7% 
(12/42) for stars with spectral types K6 and earlier, 16%±5% (10/62) for K7-M2 stars, and 
19%±4% (18/94) for M2-M6 stars. Our result is inconsistent with a disk fraction of wTTs 
peaking around K6-M2 stars and suggests that the disk fraction dependence on spectral type 
found by Lada et al. 2006, if real, is likely to be driven by the cTTs population, which would 
imply that the duration of the accretion phase is a function of spectral type. Given the very 
different spectral type distribution of our sample (there are only 7 objects with types later 
than K0 in our sample), our result does not contradict those presented by Carpenter et al. 
(2006). 

3.2.4- Fractional Disk Luminosity 

The ratio of the disk luminosity to the stellar luminosity, L disk I L*i is a measurement 
of the fraction of the star's radiation that is intercepted and re-emitted by the disk plus 
any accretion luminosity. This quantity is intimately related to the evolutionary status of 
a circumstellar disk. On the one hand, the primordial, gas rich, disks around cTTs have 
typical L DISK /Lz values > 10 — 20%. On the other hand, L DIS k/L* values for optically thin, 
gas poor, debris disks around MS stars range from 10" 3 to 10" 6 (Beichman et al. 2005). 
In order to characterize Loisk/L^ for wTTs disks, we estimate this quantity for all the 
wTTs disks in our sample except for the few disks in IC 348 for which 24 /im fluxes are not 
available. These objects are below our 24 /zm sensitivity limits, and the disk luminosities 
obtained from the IRAC excesses alone would be highly uncertain. 
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We estimate the disk luminosity according to the following procedure. We first calculate 
the IR-excess at IRAC wavelengths and at 24 /zm by subtracting from the observed fluxes 
the expected photospheric contributions predicted by the stellar models shown in Figures 
2, 3, and 4. Then, the IR-excess at 24 fim is extrapolated to longer wavelengths assuming 
the emission of a black body peaking at 24 /xm (i.e., T = 121 K) and diluted by an emis- 
sivity proportional to A -1 . Also, the IR emission of the shortest significant IRAC excess is 
extrapolated to shorter wavelengths assuming black body emission. The temperature of the 
adopted black bodies ranges from 1400 to 500 K depending on the IRAC band from which 
the IR-excess is extrapolated. Finally, we calculate the total disk luminosities by integrating 
the observed and extrapolated IR-excesses over frequency. Similarly, the stellar luminosities 
are calculated by integrating the fluxes of the stellar models shown in Figures 2, 3, and 4 
over frequency. The distribution of Ldjsk/L* for our sample of wTTs is shown in Figure 
8. The Ldjsk/L* values calculated using the same procedure for a sample of cTTs (From 
Cieza et al. 2005) and debris disks (from Chen et al. 2005) are shown for comparison. The 
Ldisk/L* values we derive for our sample should be considered lower limits because our 
assumption that the flux density peaks at 24 /im may underestimate the flux contribution 
of cool material in the outer disk. However, since none of our disks are detected at 70 /im, 
we can constrain the flux contribution of the outer disks that may be missing in our disk 
luminosity estimates. For each object, we approximate the outer disk as a diluted black 
body peaking at 70 /im (i.e., T = 41 K) with a flux density equal to 45 mJy, the estimated 
3-cr limits of the c2d observations at 70 /mi (Evans et al. 2005). For objects in Lupus and 
Ophiuchus, this approach constrains the luminosities of the outer disks to the 10 -5 -10 -4 L* 
range. For objects in IC 348, this range is 10~ 3 -10~ 4 L*. Thus, we conclude that, for objects 
with warm disks for which we derive Ldisk/L* > 10~ 3 , the flux contribution of the outer 
disks is a small fraction of the total disk luminosity. For the objects we derive Ldisk/L^ < 
10~ 3 , the contribution of an unseen outer disk could, in principle, dominate the total disk 
luminosity. However, the disk luminosities of such cold disks would still remain in the debris 
disk range. 

Figure 8 shows that only one object, Sz 96, has a Ldisk/ ' L* value characteristic of cTTs. 
This object is a M2 border line cTTs/wTTs with a Ha EW of 11 A. As shown in Figure 2, 
the SED of Sz 96 shows strong IR-excess at 2MASS wavelengths and is indistinguishable from 
that of a cTTs. This object is likely to be one of the actively accreting interlopers discussed 
in Sections 12.21 and 13.2.21 Similarly, only 4 objects have 10~ 3 > Ldisk /L* > 10 -6 , the range 
characteristic of optically thin debris disks. However, we note that the low luminosity and 
low optical depth of a disk do not necessarily warrant debris disk status, which requires the 
presence of second generation of dust produced by the collision of planetesimals in a gas poor 
environment. The lowest luminosity, and presumably optically thin, disks in our sample are 
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RXJ1622.6-2345, ROXs 36, IC348-56, and IC348-124. We model the SEDs of these disks in 
Section 14.31 If these objects are in fact debris disks, at an age of ~l-3 Myrs, they could 
be some of the youngest debris disks observed to date. However, the confirmation of debris 
disk status would require information on the grain size distribution and gas content of their 
disks. The main conclusion that can be drawn from Figure 8 is that the bulk of the wTTs 
disks have L disk I ' L* values that bridge the gap between the cTTs and the debris disk range. 
This supports the idea that wTTs disks represent an intermediate evolutionary stage linking 
primordial disks around cTTs and debris disks around MS stars. 

4. Discussion 
4.1. Comparison to recent Spitzer results 

4.1.1. Comparison to Padgett et al. 2006 

The IRAC disk fractions we found in Section 13.11 are ~3-6 times larger than those re- 
cently found for a sample of 83 relatively isolated wTTs (Padgett et al. 2006, P06 hereafter). 
P06 find that 3 of their objects show IR-excess both at IRAC and MIPS wavelengths, while 
2 objects show IRAC fluxes consistent with photospheric emission and small excesses at 24 
/im. Most of these stars in the P06 sample are members of the extended population of Li rich 
wTTs discovered by the ROSAT X-ray satellite around nearby molecular clouds and are, for 
the most part, located outside the high extinction regions mapped by the c2d project, but 
within ~6 degrees of the centers of the c2d cloud maps. The P06 observations were sensitive 
to the stellar photospheres of the entire sample at both IRAC wavelengths and 24 /im. Given 
the size of the samples involved, the discrepancy in the IRAC disk fraction of the off-cloud 
wTTs studied by P06, 4±2%, and that of the on-cloud wTTs studied in this paper, 20±3%, 
is significant at the 4-<r level. 

Follow up spectroscopic studies of ROSAT targets (e.g. Covino et al. 1997; Martin 
et al. 1999; Wichmann et al. 2000) used the strength of the Li I 6707 absorption line to 
discriminate bona-fide PMS from young active MS stars (possibly ~100 Myrs old) sharing 
the same X-ray properties. P06 selected their sample from ROSAT wTTs with Li I 6707 
absorption lines stronger than that of a Pleiades star of the same spectral type. These 
objects were suspected to be 1-10 Myr-old because the age distribution of the ROSAT 
wTTs, derived from their position in the H-R diagram, is not significantly different from 
that derived for the on-cloud TT Tauri stars (Alcala et al. 1997; see also Section 14.41 of 
this paper). This age estimate assumes that the ROSAT sources are located at the same 
distances as their adjacent molecular clouds. However, the distance and pre-main sequence 
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status of the extended population of ROSAT sources has been called into question. Feigelson 
(1996) and Briceno et al. (1997) argue that the distributed population of ROSAT sources 
consist of mostly foreground post-T Tauri and active young (~100 Myrs old) MS stars which 
are not necessarily associated with regions of current star formation (i.e. they are old enough 
to have moved far away from their birth sites) and that the presence of Li I 6707 does not 
warrant PMS status since it is highly dependent on stellar mass, convection, and angular 
momentum history. The low disk fraction of the distributed population of ROSAT sources 
studied by P06 with respect to that of the sample of on-cloud wTTs studied here supports the 
idea that the ROSAT sources found around molecular clouds represent a significantly older 
population of stars than that of the < 1-10 Myrs old cTTs and wTTs found within molecular 
clouds. The age estimate of the extended population of Li rich ROSAT sources depends on 
whether or not the P06/ROSAT targets are at the same distances as their adjacent molecular 
clouds. The luminosities of stars closer to us than the molecular clouds themselves could be 
considerably overestimated, translating into underestimated ages when the stars are placed 
in the H-R diagram. If the P06 sources are in fact associated with the molecular cloud and 
the sample has been biased toward foreground stars, the wTTS in P06 sample could be 
somewhat older than their nominal age of 1-10 Myr. It is reasonable to assume that such 
a bias toward foreground objects exists because P06 favored bright objects when selecting 
their sample in order to increase the number of wTTs that could be observed at a given 
sensitivity. If the P06 sources are not associated with their adjacent molecular clouds, then 
the only age constraints are Li I 6707 EW and X-ray brightness, and ages of the order of 
100 Myrs could not be ruled out. 



4-1.2. Comparison to Lada et al. 2006 

In a recent study of ~300 confirmed members of IC 348 (wTTs and cTTs) Lada et 
al. (2006, L06 hereafter) find an IRAC disk fraction of ~35% among IC 348 members 
classified as wTTs. The sample of wTTs studied by L06 is virtually the same sample of IC 
348 wTTs included in this paper, for which we derive an IRAC disk fraction of 23% ± 4%. 
In what follows, we attempt to account for the difference in these results. L06 use a disk 
identification criterion based on the slope, a, of a power law fit to the four IRAC bands. From 
the comparison to disk models, they identify objects with a > -1.8 as optically thick disks. 
Based on the predicted slope of an M0 star and the typical uncertainty in the power law fit, 
they identify objects with a < -2.56 as stellar photospheres and objects with intermediate 
slopes, -2.56 > a >-1.8. as "anemic disks". Cieza & Baliber (2006) suggest that many of the 
objects identified as "anemic disks" by the L06 are more consistent with stellar photosphere 
of late M stars than with circumstellar disks. Based on a large sample (> 400) of cTTs 
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and wTTs from the literature with photometric uncertainties smaller than 0.1 mag, Cieza & 
Baliber (2006) find that cTTs and wTTs disks occupy a well defined locus in the [3.6]-[8.0] 
vs [3.6]-[5.8] diagram with [3.6]-[8.0] > 0.7, while bare stellar photospheres have [3.6]-[8.0] 
< 0.5. In fact, only ~1 % of the their sample have 0.7 > [3.6]-[8.0] > 0.5. As seen in Figure 
6, faint objects in IC 348, most of which are M4-M7 stars, have very uncertain [3.6]-[8.0] 
colors. 

We suggest the possibility that the combined effect of slightly redder stellar photospheres 
and larger photometric uncertainties in late M stars with respect to brighter stars with earlier 
spectral types resulted in a large fraction of M4-M7 stellar photospheres being classified as 
"anemic" disks by L06. L06 find that the ratio of "anemic" disks to optically thick disks 
increases from ~l/5 to ~1/1 from K6-M2 to M2-M6 stars, which is consistent with the idea 
that some of the anemic disks around late type stars may be of questionable significance. 

We find that most of the wTTs that do not satisfy our disk identification criterion 
discussed in Section [37TI but are classified as "anemic" disks by L06 are objects with spectral 
types M5 or later, which tend to have very uncertain IRAC colors due to the strong IR 
background of the IC 348 cluster. These low-mass objects account for the difference in the 
disk fraction derived by us and L06 for the wTTs in IC 348. In Table 4 we list the 18 
objects classified as anemic disks by L06 that do not satisfy our disk identification criterion 
discussed in Section [37T1 We find that, for objects in Table 4, the [3.6]-[8.0] colors obtained 
by L06 are ~0.15 mag redder in the mean than our colors. However, this is a selection effect 
due to the fact that we are only considering the objects that have [3.6]-[8.0] colors from 
L06 that are red enough to be classified as anemic disks. In fact, when the entire sample 
of IC 348 wTTs is considered, our photometry gives [3.6]-[8.0] colors that are ~0.05 mag 
redder than those obtained by L06. We note that even using the colors from L06, only one 
of the objects satisfies [3.6]-[8.0] > 0.7, the disk identification criterion proposed by Cieza 
& Baliber (2006). Also, for 12 of the 18 objects listed in Table 4, L06 derive a's that are 
consistent with the expected stellar photosphere within the errors (i.e. (a+2.66)/a a < 3, 
where -2.66 is the slope of an M0 photosphere adopted by L06). Finally, two of the objects for 
which (a+2.66)/cr Q , > 3 have no 8.0 fim fluxes measured by us or by L06, which renders the 
disk identification less reliable. We conclude that an accurate identification of weak IRAC 
excesses requires careful consideration of the photometric uncertainties involved, especially 
in regions with strong extended emission such as the IC 348 cluster. 
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4.2. The diversity and evolutionary status of wTTs disks 

For late type stars, the IRAC and MIPS fluxes originate from the inner ~20 AU of the 
disk. For this reason, observations at these wavelengths are especially suited to study the 
evolution of the planet-forming region of the disk. According to models of the evolution of 
circumstellar disks, dust growth and settling take place very efficiently in the circumstellar 
disks soon after their formation (Weidenschilling 1997, Dullemond & Dominik 2005). Nu- 
merical models also predict that the growth of dust will reduce the opacity of the disk to 
the stellar radiation causing the disk to become flatter and the excesses in the mid infrared 
to decrease with time proceeding from the shortest to the longest wavelengths. There is 
mounting evidence of such phenomena in the IRAC colors of T Tauri stars and lower mass 
objects (D'Alessio et al. 2001, Furlan et al. 2005, Sicilia-Aguilar et al. 2005, Muzerolle et 
al. 2006). 

In this section, we try to explore observational consequences that could be related to 
dust grain growth and/or settling toward the mid-plane for the sample of wTTs with disks. 
To this end, we estimate the longest wavelength without significant infrared excess from the 
dereddened SEDs, A tU rn-ofr (in /-tni), and the slope of the SED longward of A tU rn-ofr, cWess 
for all the wTTs stars with disks in the sample. We compare these values for our sample of 
wTTs to those we obtain for a sample of cTTs in Chamaeleon from Cieza et al. (2005) and 
a sample of debris disks from Chen et al. (2005). Figure 9 shows the ranges of mid-IR slopes 
^excess versus the wavelength of disk emission turn off A tU rn-off- Systems with A tU rn-off < 2.0 
fim show some excess in the near-IR JHK bands and are probably actively accreting (see e.g. 
Hartmann et al. 2005 and Muzerolle et al. 2003). Objects with A turn _ g- in the IRAC range 
are likely to be purely irradiated disks with increasing degree of settling and/or larger inner 
holes with increasing A turn _ fr- Finally, those stars with A turn _ fr > 8.0 fim are likely to have 
cleared inner disks with radii of several AU (e.g. Calvet et al. 2002) or alternatively have 
dust in the inner disk which has grown to sizes large enough not to produce any detectable 
excess above the photospheric level. 

A simple disk evolution scenario based solely on dust grain growth and settling predicts 
smaller access slopes for longer A tU rn-off or larger cleared inner regions (Dullemond & Dominik 
2005, D'Alessio et al. 2005). However, Figure 9 suggests that the range of possible excess 
slopes increases with the process of inner disk clearing. The actively-accreting cTTs stars 
all cluster around access values of —1. For these objects, the onset of the IR excess occurs 
near 2.2 /xm. We find that wTTs show a richer distribution of SED morphologies than cTTs. 
Some of the wTTs show SEDs that are indistinguishable from those of cTTs, while others 
show a wide range of A tU rri— off ^excess values. Finally, the debris disks, with ages 12 and 200 
Myrs (Chen et al. 2005) show detectable excess only longward of 8.0 jum, but also present a 
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great variety of spectral slopes. 

In general, the diagram suggests an evolutionary sequence in which most actively ac- 
creting cTTs have similar near to mid-IR SEDs, dominated by optically thick emission of 
the inner disk. This similarity is probably due to the fact that, in most cases, their inner 
disks extend to the dust sublimation radius (Muzerolle et al. 2003). On the other hand, 
if disks evolve from the inside out, wTTs disks are likely to have a wider range of of inner 
disk radii and temperatures than cTTs. This diversity of inner disks, together with possible 
grain grown and grain settling to the disk mid-plane may explain the large ranges of a excess 
and A turn _ fr found. Finally, the debris disks populate the right part of the plot with large 
cleared inner disks and emission from cool optically thin dust. 



4.3. Comparison to Optically Thin Disk Models 

As discussed in Section I3T2"| four stars in our sample show SEDs and L D jsk/L* consistent 
with optically thin disks. These are ROXs 36, RXJ1622.6-2345, IC348-56 and IC348-124.In 
this section, we modeled these four objects using the the optically thin disk model developed 
by Augereau et al. (1999). 

We limited the exploration of the parameter space to the disk parameters that affect 
most the global shape of an SED, namely the minimum grain size a m j n , the peak surface 
density position r and the total dust mass Md us t (or, equivalently, the surface density at 
r ). We adopted a differential grain size distribution proportional to a -3 5 between a m i n and 
a max) with a max = 1300 jttm, a value sufficiently large to not affect the SED fitting in the 
wavelength range we consider. Following Augereau et al. (1999), the disk surface density 
S(r) is parametrized by a two power-law radial profile S(r) = T,(ro)\/2 (x~ 2ain + x~ 2a ° ut ) 1 ^' 2 
with x = r/r , and where a- in = 10 and a ou t = —3 to simulate a disk peaked around r , with 
a sharp inner edge, and a density profile decreasing smoothly with the distance from the star 
beyond r . The optical properties of the grains were calculated for astronomical silicates 
(optical constants from Weingartner & Draine (2001)), and with the Mie theory valid for 
hard spheres. The grain temperatures were obtained by assuming the dust particles are in 
thermal equilibrium with the central star. NextGen model atmosphere spectra (Hauschildt 
et al. 1999) scaled to the observed dereddened K-band magnitude, were used to model the 
stellar photospheres. 

For each star, we calculated 15000 SEDs (0.3 fim< A < 950 /im), for 75, logarithmically- 
spaced values of a min between 0.05 /im and 100 /zm, and for 200 values of r , logarithmically- 
spaced between 0.02 AU and 300 AU. For each model, the dust mass was adjusted by a least- 
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squares method, assuming purely photospheric emission (within the uncertainties, 10%) in 
the four IRAC bands, and by fitting the measured MIPS 24 /zm flux density. Models with 
flux densities in the MIPS 70 /im and 160 /xm bands, larger than the estimated 2a upper 
limits of 0.03 Jy and 0.5 Jy, respectively, were eliminated. The results are summarized in 
Table 5, and the SEDs are displayed in Figure 10. Results in Table 5 are listed for two 
different regimes of minimal grain sizes, namely, a min > 10 /im and a min < 0.5 /im. The 
first regime accounts for a scenario where all grains have grown to sizes > 10 /im and where 
smaller grains are not replenished by collisions, while the second regime accounts for a case 
where grains smaller than 0.5 jum remain in the system, either because they are collisionally 
replenished, or because pressure forces are inefficient to expel grains. However, assuming 
a gas poor environment (i.e. the dust dynamics is not controlled by the gas), dissipation 
processes acting on small grains will limit the size of the smallest grains that are likely to 
populate the disk. For ROX 36, an A2 star, the blow-out size is of the order of a few microns. 
Thus, grains significantly smaller than 0.5 /im are unbound and unlikely to be present in 
the disk. For the three M-type stars, radiation pressure is not high enough to overcome the 
gravitational force and expel grains. In this case, and provided the star mass loss rate is 
high enough, the pressure force expected from the stellar wind will set the lower limit for 
a min (Augereau & Beust 2006, Strubbe & Chiang 2006). 

In most cases, neither the position of the peak surface density r , nor the minimum 
grain size a min , can be uniquely determined with so few observational constraints, but some 
models can be eliminated. In particular, the lack of excess in the IRAC bands imposes the 
disk to be significantly dust-depleted within ~1 AU from the star. The best fits to the SEDs 
of the three low-mass stars (RXJ1622.6-2345, IC348-56 and IC348-124) are indeed obtained 
for tq values around 2-5 AU, while in the case of ROXs 36, the inner hole could be ten 
times larger. The properties of the ROXs 36 disk are in fact those of a typical j3 Pic-like 
disk: a dust mass in the 10 _4 -10 _2 M® range, a low fractional luminosity Lir/L* of about 
10~ 4 , and typical minimum collision time-scales (as calculated using the formula given by 
Backman et al. 1993) one to three orders of magnitudes smaller than the star age. Also, 
as it is the case for virtually all (3 Pic-like disks, the disk around ROXs 36 seems to be 
collisional dominated because collisions occur much faster than Poynting-Robertson (P-R) 
drag (as calculated using the formula given by Augereau & Beust 2006) for all grains larger 
than the blow-up size. The best fits to the observed excesses around the three low-mass 
wTTs (RXJ1622.6-2345, IC348-56 and IC348-124) are also consistent with optically thin 
disks observed around other young M-type stars such as AU Mic (Liu et al. 2004). These 
three objects have low luminosity ratios, L IR /L* < 2 x 10~ 2 , and dust masses between 10~ 4 
and 10~ 3 M ffi , with an upper limit of about 0.5 M®. The collision time-scales for these three 
objects are extremely small (of the order of 100 yr). 
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The preliminary models described above suggest that the disks around ROXs 36, RXJ1622.6- 
2345, IC348-56, and IC348-124 could be younger analogs of the (3 Pic and AU Mic debris 
disks, and thus some of the youngest debris disks ever observed. However, this interpre- 
tation depends on the assumption that these four wTTs disks are gas poor (as their older 
counterparts are). The timescale over which gas clearing occurs is still poorly constrained 
observationally and it is even unclear whether gas and dust are lost simultaneously as disks 
evolve from the massive optically thick to the debris disk phase. High-resolution Spitzer-IRS 
observations, such as those presented by Pascucci et al. (2007), to constrain the amount of 
gas present in these extremely young (1-3 Myrs) optically thin disks could provide crucial 
information on the gas evolution in the transition from the primordial to the debris disk 
stage. 

In Section 13.11 we found that most wTTs (~80 %) show no evidence for a disk. Since 
Spitzer observations are capable of detecting very small amounts of dust in the planet- 
forming regions of the disk (r ~0.1-10 AU), the absence of mid-IR excess imposes very 
stringent limits on the amount of dust available for planet formation around "disk-less" 
wTTs (i.e. most wTTs). In this section, we use the optically thin disk models discussed 
above to constrain the maximum amount of dust that could remain undetected within the 
first few tens of AU of the wTTs in our sample. Since the 24 /im observations are the most 
constraining, we restrict our analysis to Lupus and Ophiuchus objects that have 24 /im fluxes 
consistent with stellar photospheres and exclude the objects that remain undetected at this 
wavelength. Objects in IC 348 are also excluded because we detect the stellar photospheres 
of only a couple of them. 

Since none of our "diskless" wTTs are detected in the c2d 70 and 160 jum maps, for 
the purpose of our models, we adopt nominal upper limits of 15 and 250 mJy at 70 and 160 
/im, respectively. For each model we calculated the mass encompassed within a radius r, 
as a function of this radius. With this approach, we can estimate the maximum dust mass 
in the inner regions of the wTTs with no detectable, or marginally detected, emission in 
excess to the photospheric emission. The results are displayed in Figure [11] for the Lupus 
and Ophiuchus clouds, assuming distances to the Sun listed in Table 6. For the Lupus and 
Ophiuchus clouds, the c2d Spitzer observations constrain the observed wTTs to have less 
than a few 10 -6 M ffi of dust within 1 AU from the central star, and less than a few 10~ 4 M e 
within 10 AU. These mass upper limits, obtained for minimum grain sizes between 10 and 
100 /im, drop by about an order of magnitude when a m i n < 0.5 /im.The Lupus and Ophiuchus 
wTTs with no (or marginal) excess at 24 /im, have then inner disks that are strongly depleted, 
and only extremely cold disks with large inner holes are still theoretically possible because 
of the relatively large MIPS 70 /im and 160 /im upper limits. Such belts would resemble the 
dust rings resolved about nearby young Main Sequence stars (e.g. HD 181327, Schneider 
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et al. 2006 and ref. therein). Of course, our observations only constrain the mass of dust 
and can not rule out the presence of much larger planetesimals or planets because once dust 
grains grow into larger bodies (r ^> A), most of the solid mass never interacts with the 
radiation, and the opacity function, k v (cm 2 /gr), decreases dramatically. 



4.4. Circumstellar Disks and Stellar Ages 



In Section [3711 we studied the overall disk fraction of our sample of wTTs. In this section, 
we derive stellar ages from two different evolutionary tracks and estimate the disk frequency 
as a function of stellar age. In order to derive stellar ages from theoretical evolutionary tracks, 
it is necessary to obtain the effective temperatures and luminosities of all the targets. We 
estimate the effective temperatures directly from the spectral type of the objects according 
to the scale provided by Kenyon & Hartmann (1995). We derive the stellar luminosities 
by applying a bolometric correction (appropriate for each spectral type) to the 2MASS J- 
band magnitudes corrected for extinction and assuming the nominal cloud distance listed 
in Table 6. The J-band was chosen because the effects of extinction is less important than 
at shorter wavelengths (Aj ~ 0.26 Ay) and the emission from the disk is less prominent 
that at longer wavelengths. The bolometric corrections were taken from Hartigan, Strom & 
Strom (1994) and the J-band extinction, Aj, was calculated using Aj = 1.24 x E(Rc-Ic) — 
1.24 x ((Rc — Ic) obs ~ (Rc — Ic) o) ■ The intrinsic stellar colors, (Rc-Ic)o, come from Kenyon 
& Hartmann (1995). For objects without Rc and Ic* fluxes available, we use Aj=5.88xE(J- 
K). 



4-4- 1- Estimation of Age Uncertainties 

In order to estimate the error bars associated with the ages we derive, we first estimate 
the observational uncertainties that need to be propagated through the H-R diagram. For- 
tunately, the T e f/s and luminosities of wTTs are easier to determine than those of cTTs. 
The interplay between extinction and veiling introduces a large uncertainty in determining 
stellar luminosities of cTTs, and different results are obtained depending on the band to 
which the bolometric correction is applied (Cieza et al. 2005). In fact, the luminosities 
obtained from the J-band are systematically larger (by ~35%) than those obtained from the 
Ic band. For wTTs, Cieza et al. (2005) show that the luminosities derived from the Ic 
and J bands agree to within 5%; therefore, distance is probably the dominant uncertainty 
in calculating their luminosity. The distance uncertainties listed in Table 6 translate into 
~20-30% luminosity uncertainties. Similarly, the spectral type classification of wTTs (and 
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therefore their temperatures) is usually more accurate than that of cTTs. This is because the 
spectra of cTTs are affected by veiling, and their photospheres are highly heterogeneous in 
terms of temperature due to the presence of hot accretion columns. To estimate the effective 
temperature uncertainty of the wTTs in our sample, we adopt one spectral type subclass as 
the classification accuracy. For the reasons mentioned above, we regard the T e //'s and the 
luminosities (and therefore the ages) of wTTs as being more accurate than those derived for 
cTTs using the same procedure. 

A study of the disk fraction as a function of age can only yield meaningful results if 
the intrinsic age spread of the sample is larger than the age uncertainties attributable to 
observational errors. We verify that our sample satisfies this condition by comparing the 
spread of the ages we derive to that expected solely from observational errors. Based on an 
error budget similar to the one described above, Hartman (2001) estimated that observational 
errors will introduce an age spread of cxlog(age) = 0.18 in the logarithmic age distribution 
derived for wTTs in Taurus. In Figure 12, we plot the derived age distribution of our sample 
as calculated from two different evolutionary tracks. The mean age and the age spread are 
model dependent. According the the models presented by D'Antona & Mazzitelli (1994, 
1998, D98 hereafter) @ the logarithmic age distribution of our sample can be characterized 
as a Gaussian centered around log(age)=6.3 with <r(log age) = 0.57. Similarly, according to 
the models presented by Siess et al. (2000, S00 hereafter) □, the logarithmic age distribution 
of the sample can be characterized as a Gaussian centered around log(age)=6.6 with <r(log 
age) = 0.40. The age spreads derived are 2.2-3.2 times the value attributed by Hartmann 
(2001) to observational errors, which is also shown in Fig 13 for comparison. 



4- 4-2. Evolutionary Tracks and Stellar Ages 

In order to evaluate the degree to which ages we derive depends on the models, we 
compare the individual ages derived from the D98 and S00 models. We choose these two 
particular evolutionary tracks because they provide the appropriate mass and age range and 
are both widely used, which allows a direct comparison of our results to those from other 
papers. Figure 13a shows the ages derived for our sample using both sets of evolutionary 
tracks. The error bars for every object have been calculated by propagating into the evolu- 
tionary tracks the observational uncertainties computed as described in the previous section 
(e.g. a T e ff uncertainty equal to one spectral type subclass and a luminosity uncertainty 



3 Available at http://www.mporzio. astro. it~dantona/prems. html 
4 Available at http://www-astro.ulb.ac.be^siess/database.html 
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dominated by the distance uncertainty). Stellar ages of PMS stars are very difficult to esti- 
mate due to the large observational and model uncertainties involved (Hillenbrand & White, 
2004) and they are often taken with a high degree of (healthy) skepticism. However, even 
though the the error bars in the individual ages are large, the total age spread in the sample 
is significantly larger than the typical (observational) error bar. This is consistent with the 
analysis of Figure 12. Also, even though D98 and S00 evolutionary tracks show some sys- 
tematic differences (e.g. D98 tracks yield significantly younger ages than S00 models), the 
relative ages agree fairly well. 

Figure 13a shows that ~40% of the wTTs that are both younger than 1.5 Myrs according 
to the S00 models and younger than 0.6 Myrs according to the D98 models have circumstellar 
disks. In contrast, none of the targets that are older than 10 Myrs according to the D98 or 
S00 models has a disk. The decrease in the IRAC disk fraction with stellar age is clearly 
seen Figure 13b, where we restrict our analysis to the 198 objects in the magnitude limited 
subsample discussed in Section 13.1.21 and divide the ages we derive using the D98 and S00 
models into 4 age bins. Very similar conclusions can be drawn if our sample is combined with 
that of P06. Including the P06 samples considerably increases the statistical significance of 
the last age bin. Taken both samples together, none of the ~40 stars that are older than 
10 Myrs according to either of the tracks have an IR excess indicating the presence of a 
circumstellar disk. 



4.5. Constraint on the Timescale of Planet Building 

Figure 13 suggests that circumstellar disks, as defined by the presence of IR excesses 
at A < 10-24 /im, are very rare or nonexistent around wTTs with ages >10 Myrs. This 
timescale is very similar to that obtained by studies of the frequency of circumstellar disks 
detected in the near-IR (See Hillenbrand 2006 for a review). However, our results impose 
much stronger constraints on the time available for the formation of planets than those 
provided by previous studies. Past results based on near-IR excesses always left room for 
the possibility that stars without near-IR excess had enough material to form planets at larger 
radii not probed by near-IR wavelengths. IRAS and ISO had the appropriate wavelengths 
range to probe the planet-forming regions of the disk but lacked the sensitivity needed to 
detect all but the strongest mid- and far-IR excesses in low-mass stars at the distances of 
nearest star-forming regions. Spitzer provides, for the first time, the wavelength coverage 
and the sensitivity needed to detect small amounts of dust in the planet-forming regions 
of a statistically significant number of low-mass PMS stars. In particular, the results from 
Section 14.31 suggest that 24 /im fluxes consistent with stellar photospheres constrain the 
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amount of warm dust (T ~100 K) in the disks of our sample of wTTs to be much less 
than an Earth mass. Even though our 24 jum observations are not sensitive to the stellar 
photospheres of all the targets, taking the P06 sample and our sample together, there are 
at least ~40 wTTs with estimated ages > 10 Myrs showing photospheric fluxes at 24 /im. 
This number is likely to be a lower limit because, as discussed in Section 14.1.21 we suspect 
that the ages of many of the stars in the P06 sample have been underestimated. This seems 
to imply that after ~10 Myrs wTTs have to be in a relatively advanced stage of the planet 
formation process if they are to form planets at all. 

Since the cTTs disks older than 10 Myrs are also very infrequent, 10 Myrs seems to be a 
general upper limit for the survival of primordial disks around PMS stars. This conclusion is 
also supported by recent results from the "Formation and Evolution of Planetary Systems" 
(FEPS) project. Silverstone et al. (2006, S06 hereafter) search for IRAC excesses in a sample 
of 74 young (age < 30 Myr old) Sun-like (0.7 < M*/M Q ) stars. They divided the sample 
into two age bins, 3-10 Myrs and 10-30 Myrs. S06 find IRAC excesses for 4 of the 29 stars 
in the youngest age bin and for 1 of the 45 stars in the older age bin. The FEPS objects 
were selected based on their ages (< 30 Myrs), masses (~ 1 M ), distances (< 170 pc), and 
low infrared backgrounds, without a bias with respect to their Ha EW or IR properties. 
The five objects with IRAC excess have SEDs consistent with those of CTTS, and given the 
age uncertainties, it entirely possible that they are all younger than 10 Myrs. In fact, the 
only object with IRAC excess for which S06 adopts an age older than 10 Myrs is PDS 66. 
S06 adopted an age of 17 Myrs for PDS 66 based on the mean age of the Lower Centaurus 
Crux, but its formal age is 7-17 Myrs depending on the evolutionary track used (Mamajek et 
al. 2002). Recent sub-millimeter results extend the conclusions on the survival time of the 
material in the inner disk (r < 0.1 AU) and the planet-forming region of the disk (r ~ 1-10 
AU) to the outer disk (r ~100 AU). Andrews & Williams (2005) observed over 150 YSOs in 
Taurus and found that < 10 % of the object lacking inner disk signatures were detected at 
sub-mm wavelengths. Given the high sensitivity of their survey (3<x ~10 mJy at 850 /im) 
they conclude that dust in the inner and outer disk dissipate nearly simultaneously. 



Figures 13 also suggests that the disks of some wTTs dissipate in timescales <1 Myrfl 
Are these apparently very young disk-less objects as young as the evolutionary tracks suggest 
or is their apparent youth just a product of the large age uncertainty? Several factors 
can introduce very large errors in the age determination. For instance, the luminosity of 
foreground field stars can be grossly overestimated if the nominal cloud distances are used, 



5 This assumes that all wTTs had a disk at some point of their evolution. It is possible that some disks 
dissipate very early (t << 1 Myrs), even before the star is optically revealed and can be classified as a wTTs. 
However, accretion through a disk is considered an unavoidable step of the star formation process. 
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leading to grossly underestimated ages. Also, large errors in the spectral type classification 
might lead to large errors in the extinction and luminosities. To check for these possibilities, 
we plotted the SEDs of the objects classified as the youngest disk-less wTTs in our sample. 
Their SEDs look consistent with stellar photospheres and the overall quality of the fits 
suggests that both the spectral types and the extinction corrections are reasonably accurate. 
Still, some of these objects are ~0.5 Myrs old according to the DM98 models and ~1.0 Myrs 
old according to the S00 models, and ~15 times more luminous than main-sequence stars of 
the corresponding spectral types. 

In Section 14.41 we found that the dissipation timescale or "survival time" of wTTs cir- 
cumstellar disks ranges from less than 1 to 10 Myrs. A related timescale is the transition 
timescale from optically thick accretion disks to undetectable disks. Assuming that wTTs 
disks are the link between cTTs disks and "disk-less" wTTs, then the transition timescale, 
r, can be estimated as r= jj^f x <age>, where Ntran is the number of wTTs disks, 
Npms is the total number of PMS stars (wTTs+cTTs), and <age> is the mean age of the 
sample. Adopting the wTTs/cTTs ratio of IC 348 (~3/2) and a mean age of 3 Myrs, the 
overall IRAC disk fraction in wTTs of 20% we find in Section I3TT1 implies r ~0.4 Myrs. This 
timescale is very similar to that found for the transition timescale between an optically thick 
disk and an optically thin disk by Skrutskie et al. (1990) and Wolk & Walter (1996) based on 
the number of "transition objects" , which they define as targets without K-band excess but 
with strong IRAS excesses. The fact that the transitional timescale is significantly shorter 
than the mean disk-lifetime is inconsistent with traditional viscous evolution (Hartmann 
et al. 1998) or magnetospheric clearing models (Armitage et al. 1999), which predict a 
steady disk evolution and thus similar timescales for disk lifetimes and disk dispersal times. 
Non-steady disk evolution scenarios are required to explain the short transitional timescales 
inferred after significantly longer disk lifetimes. Such scenarios include the ultraviolet-switch 
model (Alexander et al. 2006) and the presence of gap forming planets (e.g. Quillen et al. 
2004). 

Our results from Section 4.4, constrain not only the dissipation timescale of the dust 
during the planet formation process, but also the amount of second generation dust that 
is produced during this process (Section 14.31) . Numerical models presented by Kenyon & 
Bromley (2004) predict the amount of 10 and 20 /im excess as a function of time produced 
by the formation of terrestrial planets. Detailed comparison between these kind of models 
with predictive power and the new observational constraints Spitzer is now providing could 
be highly valuable for our understanding of planet formation. 

Particularly intriguing and potentially very important objects are the ~1 Myrs old wTTs 
without any measurable IR excess (for A < 24 /im) discussed above. One possible explanation 
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for the very early dissipation of their disks is that these stars have already formed planets 
through gravitational instability, which is expected to occur at extremely young ages when 
disks are most massive (Boss 2000). Another possibility for the fast dissipation of these disks 
is the presence of close companions that could have disrupted their disks. This possibility 
can be tested with a combination of radial velocity and adaptive optics observations to search 
for companions. Initial conditions could also be responsible for the early dissipation of the 
disk, although this hypothesis is not easily testable observationally. We note that the Space 
Interferometry Mission should be able to establish whether or not these very young disk-less 
wTTs do in fact harbor planets. The presence of planets around these ~1 Myrs old stars 
would set the tightest constraints to date for the planet formation timescale. 

The properties of wTTs with disks such as their age, SED morphology, L DIS kI 'L*, etc, 
strongly suggests that wTTs disks are the link between the massive primordial disks found 
around cTTs and the debris disks observed around young MS stars. They could be arguably 
the best places to study ongoing planet formation. Before the end of its mission, Spitzer is 
likely to identify hundreds of wTTs disks in nearby star-forming regions. These objects will 
most likely be the main targets of many follow-up observations. Deep far-IR and sub-mm 
observations with Spitzer, Herschel, and Alma will allow the study of the outer regions of 
these wTTs disks and estimates of their masses. Follow-up Spitzer /IRS observations of the 
10 and 20 fim silicate features will provide important information on the evolutionary state 
of the circumstellar dust around these objects. Finally, high resolution searches for H2 and 
atomic lines, such as those presented by Hollenbach et al. (2005) and Pascucci et al. (2007), 
would be highly desirable to constrain the amount of gas available for the formation of giant 
planets in wTTs disks. 



5. Summary of results 

We present a census of circumstellar disks and report the disk frequency as a function 
of stellar age for a sample of over 230 spectroscopically identified wTTs located in the c2d 
IRAC and MIPS maps of the Ophiuchus, Lupus and Perseus Molecular Clouds. Our main 
results can be summarized as follows: 

1) In Section 13.11 we find from a magnitude limited subsample of wTTs that ~20% 
of the wTTs have noticeable IR-excesses at IRAC wavelengths indicating the presence of a 
circumstellar disk. 

2) The disks frequencies we find in the 3 clouds we consider are ~3-6 times larger than 
that recently found by P06 for a sample of 83 relatively isolated wTTs projected outside the 
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boundaries of nearby molecular clouds. This discrepancy in the disk fractions supports the 
idea that samples of wTTs (nominally 1-10 Myrs old) located a few degrees away from their 
parent molecular clouds represent an older population of stars. The disk fractions we find 
are more consistent with those obtained in recent Spitzer studies of wTTs in young clusters 
such as IC 348 and Tr 37. However, in Section 14.1.21 we suggest that Lada and colleagues 
might have overestimated the disk fraction of wTTs in IC 348 by classifying as "anemic 
disks" some disk-less M4-M7 stars with large photometric uncertainties. 

3) In Section 13.2.21 we find that the disk fraction of wTTs is a smooth function of Ha 
EW. In Section f3. 2. 41 we show that the fractional disk luminosities of wTTs disks bridge the 
gap between the cTTs and the debris disk range. 

4) In Section 14.3} we estimate mass upper limits of dust within the inner 10 AU of 
10 _4 M e for the objects in our sample with 24 jum fluxes consistent with stellar photospheres. 

5) In Section I4.4[ we place our sample of wTTs in the H-R diagram and find that the 
stars with excesses are among the younger part of the age distribution. However, we also find 
that up to ~50% of the apparently youngest stars in the wTTs sample show no evidence of IR 
excess. This suggests that the circumstellar disks of a sizable fraction of pre-main-sequence 
stars dissipate before the stars reach an age of ~1 Myr. 

6) Also in Section I4.4[ we find that none of the stars in our sample apparently older 
than ~10 Myrs have detectable circumstellar disks. Since Spitzer observations probe planet- 
forming regions of the disk (r ~0.1-10 AU) and are capable of detecting IR excesses produced 
by very small amounts of dust, our results impose stronger constraints on the time available 
for the formation of planets than those provided by previous studies based on detections of 
disks in the near-IR. 

7) Finally, in Section 14.51 we estimate a transition timescale of ~0.4 Myrs between 
optically thick accretion disks and disks that are undetectable shortward of ~10 fim, in 
good agreement with previous results. 
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Table 1. Spitzer data 



ID 


R.A 


Dec 


F3.6 


Error 3 . 6 


F4.5 


Error 4 . 5 


F5.8 


Error 5 . 8 


F8.0 


Errorg. F 2 4 Error 2 4 




(J2000.0) 


(J2000.0) 










(mJy) 








IC348-1 


55.8837 


32.1048 


1.06e+01 


1.51c-01 


7.09c+00 


8.78e-02 


4.81e+00 


5.52e-02 


2.79c+01 


4.68e-02 


IC348-2 


55.8903 


32.0293 


4.96c+00 


7.22e-02 


3.42c+01 


4.89e-02 


2.52c+00 


3.64e-02 


1.42c+01 


2.89e-02 


IC348-3 


55.9526 


32.2308 


3.23e+00 


4.64e-02 


2.39c+00 


2.82e-02 


1.65c+00 


2.96e-02 


9.38e-01 


3.05e-02 


IC348-4 


55.9532 


32.1259 


9.01c+00 


1.75c-01 


6.51c+00 


8.82e-02 


4.13c+00 


5.04e-02 


2.34e+00 


4.26e-02 


IC348-5 


55.9558 


32.1778 


7.47e+00 


1.13c-01 


5.55e+00 


6.76e-02 


3.57e+00 


4.47e-02 


2.13e+00 


4.04e-02 



Note. — [The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample to illustrate 
its content.] 



Table 2. non- Spitzer data 



ID 


R.A 


Dec 


SpT 


H Q 


V Errory R c 


Error Rc 


lc 


Error/ c 


J 


H 






(J2000.0) 


(J2000.0) 




(A) 






(mag) 










IC348-1 


55.8837 


32.1048 


M0.75 


1.0 


14.81 


0.03 


13.73 


0.03 


12.18 


11.38 


11.13 


IC348-2 


55.8903 


32.0293 


M5 


5.0 


17.93 


0.03 


15.91 


0.03 


13.36 


12.57 


12.22 


IC348-3 


55.9526 


32.2308 


M5 


6 


17.92 


0.03 


16.05 


0.03 


13.59 


12.98 


12.64 


IC348-4 


55.9532 


32.1259 


Ml. 5 





15.44 


0.03 


14.16 


0.03 


12.44 


11.59 


11.34 


IC348-5 


55.9558 


32.1778 


M3.5 


11 


16.36 


0.03 


14.76 


0.03 


12.58 


11.80 


11.54 



Note. — [The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a 
sample to illustrate its content.] 
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Table 3. Adopted Extinction Relations 





A 


Ay/Ay} 


Band 


w 




V 


0.55 


1.00 


Rc 


0.65 


0.79 


lc 


0.80 


0.58 


J 


1.25 


0.26 


H 


1.66 


0.15 


K 5 


2.2 


0.09 


IRAC-1 


3.6 


0.04 


IRAC-2 


4.5 


0.03 


IRAC-3 


5.8 


0.02 


IRAC-4 


8.0 


0.01 


MIPS-1 


24 


0.002 


lr The 


extinction re- 


lations 


for the 


optical 



and 2MASS wavelengths 
come the Asiago database 
of photometric systems 

( http: / / ulisse.pd.astro.it / Astro / ADPS/ enter.html 
; Fiorucci & Munari 2003), 
while those for the Spitzer 
bands come from Huard et 
al. (2007, in prep. 



Table 4. Objects classified as anemic disks by L06 that do not satisfy our disk identification criterion 



Ra 


Dec 


ID 


SpT 


a 


(T _ 

u OL 


(a+2 66) /a„ 


[3 6148 01 


13 61-18 01 


A113 61-18 011 


(J2000.0) 


(J2000.0) 


L06 




L06 


L06 


L06 


L06 


C2D 


L06-C2D 


55.9526 


32.2308 


261 


M5 


-2.560 


0.048 


2.0 


0.29 


0.26 


0.03 


55.9742 


32.1251 


254 


M4.25 


-2.549 


0.081 


1.3 


0.30 


0.29 


0.01 


56.0185 


32.0817 


303 


M5.75 


-2.197 


0.125 


3.7 


0.61 


— 


— 


56.0740 


32.0799 


169 


M5.25 


-2.440 


0.099 


2.2 


0.41 


0.36 


0.05 


56.0816 


32.0403 


322 


M4.25 


-2.426 


0.046 


5.0 


— 


— 


— 


56.0971 


32.0318 


1684 


M5.75 


-1.811 


0.500 


1.6 


— 


— 


— 


56.1203 


32.0730 


385 


M5.75 


-2.183 


0.350 


1.3 


0.71 


— 


— 


56.1309 


32.1915 


226 


M5.25 


-2.233 


0.265 


1.6 


0.60 


0.38 


0.22 


56.1358 


32.1451 


33 


M2.5 


-2.483 


0.157 


1.1 


0.32 


0.16 


0.16 


56.1365 


32.1544 


88 


M3.25 


-2.355 


0.081 


3.7 


0.46 


0.48 


-0.02 


56.1460 


32.1269 


8024 


K7 


-2.213 


0.218 


2.0 


0.62 


0.25 


0.37 


56.1590 


32.1727 


353 


M6 


-2.251 


0.178 


2.2 


0.61 


0.16 


0.45 


56.1794 


32.1709 


217 


M5 


-2.312 


0.043 


8.0 


0.52 


0.43 


0.09 


56.1822 


32.1800 


360 


M4.75 


-2.303 


0.064 


5.5 








56.1861 


32.1251 


218 


M5 


-2.294 


0.049 


7.4 


0.54 


-0.04 


0.58 


56.1902 


32.1864 


413 


M4 


-2.398 


0.500 


0.5 




0.31 




56.2035 


32.2228 


178 


M2.75 


-2.520 


0.080 


1.7 


0.31 


0.23 


0.08 


56.2527 


32.1387 


344 


M5 


-2.376 


0.229 


1.2 


0.49 


0.45 


0.04 



Table 5. Optically thin disk properties, for two different regimes of minimal grain sizes, a 



Star 



{&min < 0.5 /an) 
r (AU) M dus4 (10- 3 M e ) 



t co «(yr) ro(AU) 



, > 10 Aim) 

(10- 3 M e ) 



(yr) L /fl /L,xl0 3 



ROXs 36 


oo+lOO 


6-7±^ 5 


77nnn +110000 

' ' uuu -74000 


i 7+2.8 
1 ' -10 


c q+5.9 
°-°-4.9 


74nn+ 1900 


o.n±g:SI 


RXJ1622. 6-2345 


4.0t& 


0.131^3 


q 1n +2000 
oiu -300 


l-9±?i 


H qc+85 


130±I? 


o.eot^ 


IC348-56 


r o+59 
°- z -4.4 


0.531™ 


lqo +1300 
lyu -190 


2-5±?;S 


i c + 180 


86^ 


1 4+ 9 - 7 


IC348-124 


4 1+ 27 - 


0.74^ 7 6 ° 9 


65l 3 f 


1 s+ 2 - 4 

1 -°-0.94 




24±« 


q 1 +8.6 

■^■■--o.es 



Table 6. Adopted Distances 



Cloud 


Distance 


Reference 




(PC) 




Ophiuchus 


125±20 


de Geus et al. (1989) 


Lupus I 


150±20 


Comeron et al. (2006), in prep. 


Lupus III 


200±20 


Comeron et al. (2006), in prep. 


IC348 


320±30 


Herbig (1998) 
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Fig. la and lb.— [3.6] vs EX([3.6]-[24]) (left) and [3.6] vs EX([3.6]-[8.0]) (right) diagrams 
for our sample of wTTs used for disk identification. See text for definitions. The 3 and 5-a 
dispersion of the stellar photospheres are shown as dashed and dotted lines, respectively. 
Most of the disks are detected at 24 fxm. 
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Fig. 2. — SEDs of wTTs disks in Lupus. The open squares represent the observed optical, 
2MASS, IRAC and MIPS-24 /im fluxes, while the dots correspond to the extinction cor- 
rected values. Model photospheres corresponding to published spectral types are shown for 
comparison. 
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Fig. 3.- 



SEDs of wTTs disks in Ophiuchus. 
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Fig. 4.— SEDs of wTTs disks in IC 348. 
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Fig. 4. — Continued. 
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Fig. 4. — Continued. 
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Fig. 5a and 5b. — The Figure on the left shows the [3.6]-[24] vs. [3.6]-[8.0] colors of our 
sample of wTTs stars. Based on this diagram, we identify three different groups: (1) stellar 
photosphere with [3.6]-[24] < 0.7 and [3.6]-[8.0] < 0.4, (2) Objects with [3.6]-[24] > 0.7 and 
[3.6]-[8.0] < 0.4 which show significant 24 /im excess but no evidence for 8.0 //m excess, and 
(3) objects with [3.6]-[24] > 0.7 and [3.6]-[8.0] > 0.4 which show evidence for both IRAC 
and MIPS excesses. Objects in the second group are likely to have optically thin disks (see 
Section 13.2.4)) . The Figure on the right combines our sample of wTTs with a sample of 
cTTs from Hartmann et al. (2005) and Lada et al. (2006). 
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Fig. 6a and 6b. — The Figure in the left shows the [3.6]-[4.5] vs. [5.8]-[8.0] colors of our 
sample of wTTs. Faint IC 348 members ([8.0] > 10.8 mag, 3 mJy) tend to have more 
uncertain colors than the rest of the sample and are shown as smaller open boxes. The 
dotted lines represent the approximate boundaries of the color of the stellar photospheres. 
Stars in the upper right corner of the diagram have both, 4.5 and 8.0 /xm, excesses. Stars in 
the lower right corner of the diagram are stars with 8.0 fim excess but no 4.5 fim excess. The 
Figure in the right combines our sample of wTTs with cTTs from Hartmann et al. (2005) 
and Lada et al. (2006). 
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Fig. 7a and 7b. — Histogram of the Ha equivalent width for stars with and without a disk 
(a), and the disk fraction of wTTs vs. Ha equivalent width (b). The disk fraction of wTTs 
seems to be a smooth function of Ha equivalent width. 
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Fig. 8. — The fractional disk luminosities, L DIS k/L*, derived for our sample of wTTs disks. 
The L DIS k/L* values calculated using the same procedure for a sample of cTTs (From 
Cieza et al. 2005) and debris disks (from Chen et al. 2005) are shown for comparison. The 
Ldisk/L* values of of wTTs disks fill the gap between the ranges observed for typical cTTs 
and debris disks, which are shown for comparison. 
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Fig. 9. — Distribution of excess slopes a ex cess vs. the wavelength at which the infrared 
excess begins A tU rn-off for the sample of wTTs (solid dots), a sample of cTTs in Chamaeleon 
from Cieza et al. (2005), the median SED of cTTs in Taurus from D'Alessio et al. (1999) in 
asterisks (marked as D99), and a sample of Debris disks from Chen et al. (2005) in diamonds. 
The diagram shows a much larger spread in inner disk morphologies of wTTs with respect 
of those of cTTs. 
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Fig. 10. — wTTs with 24 //m excess consistent with optically thin disks. On each plot, the 
darkest regions correspond to the most likely fits to the SEDs. The dashed line shows the 
thermal emission for the best-fit model, while the dotted line corresponds to the total disk 
emission (i.e. including scattered light emission.) 
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Fig. 11. — Maximum encompassed dust mass as a function of the distance from the star 
for the Lupus and Ophiuchus clouds (respectively left and right panels). The red area 
corresponds to mass upper limits when minimum grain sizes a m i n between 0.05 /im and 
0.5 /im are considered, while the blue area corresponds to 10 /im< a min < 100 /im. 




Fig. 12a and 12b. — The distributions of stellar ages obtained for our sample using the 
D98 and S00 evolutionary tracks (left and right panels, respectively). The observed age 
distributions (solid lines) are significantly wider than what is expected from the propagation 
of observational errors into derived stellar ages (dotted lines). 




Fig. 13a and 13b. — The stellar ages derived for our sample of wTTs using two different sets 
of evolutionary tracks (D98 and S00). The error bars have been calculated adopting a T e ff 
uncertainty equal to one spectral type subclass and a luminosity error calculated from the 
uncertainty in the distance (a). A clear decrease in the disk fraction is seen with increasing 
age. ~40% of the targets that are both younger than 1 Myrs according to D98 tracks and 
younger than 2 Myrs according to S00 tracks have disks. None of the stars that are older 
than ~10 Myrs according to either of the models have disks (b). 



